Gravitational waves from a dust disk around a Schwarzschild black hole 
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We investigate gravitational waves from a dust disk around a Schwarzschild black hole to focus 
on whether we can extract any of its physical properties from a direct detection of gravitational 
waves. We adopt a black hole perturbation approach in a time domain, which is a satisfactory 
approximation to illustrate a dust disk in a supermassive black hole. We find that we can determine 
the radius of the disk by using the power spectrum of gravitational waves and that our method 
to extract the radius works for a disk of arbitrary density distribution. Therefore we believe a 
possibility exists for determining the radius of the disk from a direct observation of gravitational 
waves detected by the Laser Interferometer Space Antenna. 

PACS numbers; 04.30.Db, 04.25.Nx, 04.70.-s 



I. INTRODUCTION 

Gravitational wave projects have been launched in 
many countries to explore a new world in gravitational 
physics, in astrophysics, and in cosmology. Nowadays, 
several ground-based laser interferometric gravitational 
wave detectors such as TAMA300, LIGO, GEO600, and 
VIRGO have been seeking for a direct detection of grav- 
itational waves (See Ref. 1] for the current status of 
the worldwide network of gravitational wave detectors). 
The sensitive frequency for the ground-based interferom- 
eters is within the range of 10 - 1000 Hz. The repre- 
sentative sources of gravitational waves are the inspiral 
of binary systems composed of neutron stars and black 
holes, spinning neutron stars, and low-mass x-ray bina- 
ries 12 ■ There will be also a gravitational wave detector 
in space. Laser Interferometer Space Antenna (LISA). 
The sensitive frequency band in LISA is in the range of 
10"'' — 10""' Hz, and one of the representative gravita- 
tional wave sources this detector is studying is supermas- 
sive black holes whose masses are M ~ 10^ — 1O^M0. 
There is increasing evidence that a supermassive black 
hole exists at the center of every galaxy 01 . 

There is abundant observational evidence for dense gas 
in galactic nuclei [e.g.,0, Q. For example, our galactic 
nuclei contains a 4 x IQ^Mq supermassive black hole sur- 
rounded by an ~ lO^M© molecular gas torus [e.g.,0|. 
Also the disruption of a compact star by a supermassive 
black hole spreads the gas in black hole spacetime 0| ■ 
From a theoretical point of view, the tidal disruption of 
a star by a supermassive black hole has been studied in 
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Newtonian gravity Je.g ., 0, 0, 0, ^| and in post- 
Newtonian gravity [l2|, where the black hole is regarded 
as a point mass. There are also several studies of tidal 
disruption in a Schwarzschild black hole [3 and in 
a Kerr black hole [Tfil ItR ITslIT^. In all cases, the frag- 
ments of the disrupted star are either swallowed by the 
supermassive black hole or following a highly eccentric 
orbits. There is also an indication, based on an observa- 
tion, that a disk can form around a supermassive black 
hole 113. 

Black hole perturbation approach is one of the sat- 
isfactory tools to illustrate gravitational waves from a 
compact star around a supermassive black hole. Many 
works have been published in p ost-Newtonian expansion 
|2ll |. in numerical analysis |22l | to compute gravitational 
waves from a test particle in a Kerr black hole. How- 
ever all of the above works have been investigated in a 
Fourier domain, which means that the orbital time of 
the particle is infinite. The basic equations of gravita- 
tional waves are ordinary differential ones, but we have 
to treat "infinite" orbital time by restricting our orbit 
to a special one. There is another way of computing 
linearized Einstein equations: a time domain approach, 
which directly solves the basic equations in partially dif- 
ferential ones. The main advantages of this time domain 
approach are that we can impose an initial condition as 
an appropriate distribution of the compact object and 
that we can treat an arbitrary orbit of the particle in gen- 
eral. There are many studies in the field of time domain 
computation of the head-on collision of two black holes 
perturbatively ,2^, 2As -25J, of black hole perturbation in 
Schwarzschild spacetin ie l26l l2^ l28l I29I Isof and in Kerr 
spacetim e Isll Is^. Is^. l34j . of perturbation in spherical 
stars l36l l37l l38l l39| . of perturbation in nonrotating 
objects and of perturbation in slowly rotating stars 

Our purpose in this paper is to investigate whether we 
can extract any physical property of the disk around a su- 
permassive black hole from gravitational waves. Several 



2 



studies have investigated extracting physical properties 
of a tidally disrupted star in supermassive black holes 
mis liH E Eal • When the disk is formed from tidal 
disruption of a star by a supermassive black hole, test 
particle approximation to illustrate the fragments is a 
satisfactory tool since the pressure gradient is no longer 
dominant l4^. Saijo and Nakamura ^3iE3 investigated 
a tidally disrupted star falling into a black hole and found 
that we can extract the size of the star from the energy 
spectrum of gravitational waves. Although they inves- 
tigated gravitational radiation in a Fourier domain, in- 
appropriate description of the pre-disruption of the star, 
their Fourier domain approach does almost illustrate the 
picture of post-disruption appropriately since most grav- 
itational waves are radiated at a quasi-normal ringing 
phase. Here we mainly focus on gravitational waves from 
a disrupted star around a supermassive black hole in a 
time domain to learn whether we can extract any phys- 
ical property of a dust disk. We have to deal with two 
types of test particle orbit in our dust disk calculation 
due to the angular momentum depletion. Particles both 
plunging into and orbiting a Schwarzschild black in the 
same calculation. To perform the above situation from 
the computational point of view, time domain approach 
is quite easy to handle rather than Fourier domain one. 

This paper is organized as follows. In Sec.|nl we briefly 
describe the basic equations of gravitational waves from 
a test particle. In Sec. IllII we demonstrate gravitational 
waves from a test particle in our developed code, com- 
paring with well-known results in this field. In Sec. IIVI 
we show gravitational waves from a dust disk around a 
Schwarzschild black hole. We summarize our paper in 
Sec.0 Throughout this paper, we use geometrized units 
G = c = 1, where G and c denote the gravitational con- 
stant and the speed of light, respectively, and the metric 
sign as 



II. BLACK HOLE PERTURBATION 
APPROACH IN SCHWARZSCHILD SPACETIME 

Here we describe the basic equations of gravitational 
waves in order to compute gravitational waves from a 
test particle around a Schwarzschild black hole. 

First, we describe the motion of a test particle with 
a rest mass /i around a Schwarzschild black hole with a 
gravitational mass M . In this paper, we only consider 
the motion of a test particle in the equatorial plane. The 
motion is specified by giving the two parameters E-p and 
Zp in units of fx as 



Ep/n, Lp = Lp/^, 



(2.1) 



where Ep and Lp are the energy and the total angular 
momentum of a test particle observed at infinity, respec- 
tively. With these parameters, the equation of motion of 



a test particle is described as 
df' ~ Ep'K r 



M ^ L 



r 



dcf) _ I f 2M\ Lp 
dt ~ r'^K r ) 



(2.2) 
(2.3) 



where = r + 2Af ln(r/2M — 1) is the tortoise coor- 
dinate. We also introduce the effective potential Vp by 
integrating the radial motion of a test particle (Eq. |2.2) ') 
as 



dr.^, 
~dt 



En 



(2.4) 




The potential Vp has a local minimum when Lp > 2^/3 AI 



at the radius 



Lp2 - 12M2 



(2.6) 



where the particle has a stable circular orbit. 

Next, we describe the basic equations of gravitational 
waves. These equations are obtained by linearizing 
the Einstein equations in Schwarzschild spacetime called 
Regge- Wheeler equations for odd parity Zerilli 
equations for even parity Note that we can only 

treat odd and even parity separately in the nature of 
spherical symmetric background spacetime. 

The radial wave function of gravitational waves obeys 
Regge- Wheeler- Zerilli equations as 



I X/, ('^^ (odd/oven) 

dt^ drl ' 

_ o(odd/cvon). n 



(odd/cvcn) 
Im 



(2.7) 



1 .T-(odd/cvcn) ^(odd/cvcn) T^(odd/cvcn) . ,i i. i 

where v]^ , 5;,^ j is the radial 

wave function, source term, and effective potential of 
gravitational waves, respectively. The effective potentials 
for odd and even parities are given by 



(2., 



_ 2/ (A2(A + l)r3 + 3Af A^r^ -t- QM^Ar + MI^) 

(2.9) 



where / = 1 - 2M/r, A = A -I- UL/r, and \={l + 2){l- 
l)/2. The source terms are constructed from the energy 
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momentum tensor of a test particle as 



It is written as [Eq. (3.1) of Ref. |28 



S. 



(odd) _ 
Im 



Snrf 



(even) 
Im 



v/2A(A + l) 



-r [J^lm) V/m 

dr r 



(2.10) 



(A + 1)A 
+ r/2(A-/)A„„ + A/2r/2Gi™ 

- i^[A(A - ly + (4A - 9)Mr + 15M']aI°2 



167rr/2 _ 167rr/ 



AVXTT' 



V2A(A + 1) 



(2.11) 



where 



,S{r-R) 



dR 
'dt 



Sir - R) 
{r - 2M)2 



M7 



1 dR5{r-R) dYhni^) 



dt 



\ + l dt r - 2M 
2v/2A(A + 1) \^ rfi y 

^.(.-i^)(l)%™(^7), 



G, 



I 1 6{r- R)dRd<^>dYi^{n) 
A(A + 1) r-2M 'dt'dt 80 ' 



(2.12) 
(2.13) 
(2.14) 
(2.15) 

(2.16) 
(2.17) 
(2.18) 



Note that A^^, Ai,n, Bi,n, Di^, Fim, G^n, and are 
the coefficients expanded in tensor harmonics j4il |. Note 
also that 7 = dt/dr = Ep/ (1— 2Af/r), where r is a proper 

time of the particle. X^{t) ~ (^R-it), f^(i)^ is the location 

of a particle, where Cl = ^{t)j . Also a quantity A 

represents the complex conjugate of a quantity A. The 
functions Xim and Wim are the tensorial part of tensor 
harmonics defined by 



Xirrt = 2d4,{de - cot^)!^^, 



dl 



Yi 



sm 



Im • 



(2.19) 
(2.20) 



Finally, we mention the relationship between the radial 
wave function and the perturbed metric in the wave zone. 



h+-ih^ = —J2 



1 {1 + 2)1 



-2Yi^{e,ip) 



2rf;^J [1-2)1 
x(<™^^-2^|*d^'4r)) ,(2.21) 

where spin- weighted spherical harmonics -2^im(^, V^) is 
defined by 



-2Ylm{0,(p) 



2V^(A + 1) 



smt 



Xi^ . (2.22) 



The energy flux [dE^^/dt) of gravitational waves at the 
wave zone is also given by [Eq. (B.26) of Ref. JSJ] 



dE, 



gw 



1 ^ {l + 2)\ 
dt 647r^(?-2)! 



9 T 



dt 



(even) 
lm 



(odd) 



(2.23) 



III. GRAVITATIONAL WAVES FROM A TEST 
PARTICLE IN SCHWARZSCHILD SPACETIME 



Here we compute gravitational waves from a test par- 
ticle in Schwarzschild spacetime. We use a time domain 
approach to solve Eq. (|2.7|) . In order to test our newly 
developed code, we compare our results in two cases 
with well-known semi-analytical results. One is gravi- 
tational waves from a particle in a circular orbit around 
a Schwarzschild black hole, which has been studied up to 
the 5.5 post-Newtonian order The other is gravita- 
tional waves from a particle with a radial in-fall into a 
Schwarzschild black hole in a Fourier domain • 

We set the initial condition of the radial wave function 

as 



(odd/c 



= 0, 



^odd/cvcn) 
Irn 



t=0 



dt 



0, (3.1) 



with the boundary condition that there is no incoming 
wave to our system. Namely, we impose the outgoing 
wave boundary condition for the radial wave function at 
horizon and at a spatial infinity as 



d_ _d_ 

dt 9r, 
d___d_ 
dt 



(odd/cvcn) 



at n = (3.2) 



(odd/c 
im 



= at 



\ (3.3) 



where r™ and r°^^ are the inner and outer edge of the 
computational grid. In practice Eqs. H3.2|) and H3.3|l are 
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written as 



Circular orbits 



(odd/even) 
Im 



(t + M, 



^(odd/oven) 
Im 



At 



(odd/cvcn) 



+ 



At 
Arl 



(odd/cvcn) 



(t + At,rr) = (l-^)* 



- (odd/c 



' + An), 
(3.4) 

^(t,rr) 



At 
A^ 



(odd/c 
Im 



\t,rT'-Ar,), 
(3.5) 



We compute the radial wave function from a test par- 
ticle in a circular orbit at Rc = lOM in Fig.^ where -Rc 
is the radius of the orbit. We set the standard deviation 
parameter as a = O.IM and set the location of the ob- 
server at i?obs = 2000M. We only show the I — m — 2 
mode, which is the dominant one for the observer in the 
equatorial plane. We find a burst in the early stage of 
the evolution, which comes from the inappropriate ini- 
tial condition we imposed (Eq. j3.1j 'l. Since the source 
term of Eq. H2.7|l is unaffected by our initial condition, we 
simply ignore the unphysical burst in the wave function. 



where At and Ar* are the step-size for time and space 
respectively. Note that we apply the first order finite dif- 
ferencing scheme in both space and time to the boundary 
conditions. 

Although the initial condition we choose (Eq. |3.1j ') 
is somewhat artificial, the particle motion is unaffected 
by gravitational radiation in our test particle approxima- 
tion, and hence, after a certain evolution time, the radial 
wave function should converge into an appropriate solu- 
tion. Therefore, we simply ignore the early stage of the 
evolution throughout this paper. 

We should carefully deal with the delta function 
6{r — R), which appears in the source term [Eq. (|2.7I) ]. 
One treatment is to approximate the delta function to 
a step function that only covers one grid cell in the null 
coordinates 0, The other is to adopt, as we do 

here, an approximate function to describe a delta func- 
tion, which appears in the source term, as [s^ 



6{r -R) = 



■ exp 



(r - i?)2 



(3.6) 



where cr is a standard deviation parameter. Note that 
Eq. (|3.6|) becomes exactly a delta function in the limit of 
cr -> 0. 

We solve Eq. (|2.7|) in second order accuracy in space 
and time with the scheme based on two times iterated 
Crank-Nicholson method [Ull^l. We set the computa- 
tional grid as -lOOOM < r* < i?obs + 2000M, where i?obs 
is the location of the observer in the equatorial plane. We 
also set the step-size for space and time as Ar, = l.OM 
and At = 0.5 x Ar*, which satisfies the Courant condi- 
tion. Although we vary the location of the inner bound- 
ary < — lOOOM, the difference of gravitational wave 
amplitude is within the round-off error. 

We set the standard deviation parameter cr as ct = Ar^, . 
When cr is less than Ar*, it is not certain to specify the 
location of a test particle due to the less spatial grid 
resolution, and therefore we discard such region. More 
concrete examinations of the dependence of cr on the cir- 
cular and free fall orbits are discussed in the next two 
subsections. 
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FIG. 1: Radial wave function of a test particle in a circular 
orbit at Rc = lOM. We only show the I = m = 2 mode, 
which is the dominant mode according to the observer in the 
equatorial plane at -Robs = 2000M. We set the standard de- 
viation parameter as cr = O.lAf. Solid and dash-dotted lines 
represent real and imaginary parts of the radial wave function, 
respectively. 



We vary the location of the observer i?obs to investi- 
gate the radius which is numerically regarded as infinity 
in Fig. 121 We only compute I — m — 2 mode, and vary 
the standard deviation parameter cr. It is clearly shown 
in Fig. 12 that the amplitude becomes almost constant 
when i?obs > 2000M, irrespective of the standard devi- 
ation parameter a. Hereafter, we set the location of the 
observer at i?obs > 2000M. 



We compare our numerical results with the ones of 
5.5 post-Newtonian expansion |50l | at two different radii, 
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0.250 



0.245 



0.240 



0.235 



a = l.OM 
cr = 0.5M 
cr = 0.25M 
a = 0. IM 



1000 2000 3000 4000 5000 
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TABLE I: Comparison of our numerical results with the re- 
sults of 5.5 post-Newtonian expansion (Eq. 13.71 '). 



Re" 


a" 


dEg^/dt [i?/M^Y 


relative error [%] 


7M 


LOOM 


4.15 X 10"^ 


5.87 


7M 


0.50M 


4.04 X lO"'^ 


3.06 


7M 


0.25M 


4.01 X lO"'^ 


2.30 


7M 


O.IOM 


4.00 X 10"^ 


2.04 


lOM 


LOOM 


6.35 X lO"'^ 


3.59 


lOM 


0.50M 


6.20 X 10"^ 


1.14 


lOM 


0.25M 


6.16 X lO"'^ 


0.489 


lOM 


O.IOM 


6.15 X 10"^ 


0.326 



"i?c: Radius of the orbit 

^a: Standard deviation parameter 

'^dEg„/dt: Energy flux of gravitational waves 



FIG. 2: The amplitude of the radial wave function as a func- 
tion of the location of the observer iiobs. We only show the 
I = m = 2 mode. Solid, dash-dotted, dashed, dotted lines 
denote the standard deviation parameters a — l.OM, 0.5M, 

0. 25M, O.IM, respectively. 

1. e., Rc — 7M and lOM. The energy flux of gravita- 
tional waves {dEg^ /dtY^^'f from the analytical calcula- 
tion (Eq. [3.1] of Ref. IH^) is that 

/ dEgw \ ^™'> _ ( 3.92 X IQ-^ifi/Mf for Rc = 7M. 
\ dt ) ~ { 6.13 X 10-%i/M f for Rc = lOM. 

Our energy flux for four different standard deviation pa- 
rameters (7 at i?obs = 2000M is summarized in Table 
We sum up the index of the spherical harmonics / from 2 
to 6. From Table U we find that the deviation from the 
results of 5.5 post-Newtonian expansion becomes small 
when the standard deviation parameter a becomes small. 
It is certainly true that this difference disappears within 
the limit of u 0. We therefore conclude that for a 
circular orbit our code has good accuracy for a ^ 0.5M 
within the relative error of few percent in the energy flux. 

We also mention the convergence order of our code. 
Since we use finite differencing scheme and approximate 
the delta function, there are two parameters a and Ar, 
for investigating the convergence order. 

First we fix a as 0.5M and varies Ar^, between O.IM 
and 0.5M. The energy flux (dEg^/dt) at Rc = lOM are 
6.17X 10-5[AiVA/2] for Ar, = O.IM, 6.17 x 10"5[/iVM2] 
for Ar* — 0.2M (relative error from the energy flux of 
Ar* = O.IM is 0.07%), 6.19 x 10-^[fi^/M^] for An = 
0.4M (relative error 0.35%), and 6.20 x lO-^i/x^/M^] for 
Ar* = 0.5M (relative error 0.56%). Therefore our code 
roughly has second order convergence in space, which 
is consistent with our second order flnite differencing 
scheme. Next we flx the relationship a = Ar* and varies 



a between O.IM and l.OM. Roughly speaking our code 
has flrst order convergence in a, since the relative error 
in Table n decreases one half as a becomes one haff. 

B. Free fall orbits 

Next we compute gravitational waves from a test par- 
ticle falling into a Schwarzschild black hole. We show the 
radial wave function oil = m = 2 mode at i?obs = 7000Af 
in Fig. 121 Instead of following the geodesies of a parti- 
cle from infinity, we put an in-fall velocity at r = rjnt of 
V = dr^/dt = — -\/2M/rint derived from Eqs. H2.4|) and 
(|2.5() . We also vary the initial radius of the test par- 
ticle Tint in the range of 200M < nnt < 3000M, and 
find that the total amount of radiated energy is in the 
relative error range of 0.001%. The radiated energy is 
determined by the integration of Eq. (|2.23() in the range 
of 1600M <t-r< 3000M in Fig.EI We set the spatial 
grid size Ar* — o-, as we mentioned before. We show our 
convergence test of the ringing tail, varying Ar^ in Ta- 
ble ^1 We find that our computation has a convergence 
when we take a < 0.25Af. The resonance frequency 
agrees with the early results of Sasaki and Nakamura 
[54l | (Muj ~ 0.30), which corresponds to the frequency of 
the quasinormal mode of a Schwarzschild black hole. 



TABLE II: Resonance frequencies of the quasi-normal ringing 
for some values of Ar* and a. 

Ar, a Muj 
O.IOM O.IOM 0.31 
0.25M 0.25M 0.30 
0.50M 0.50M 0.29 
LOOM LOOM 0.20 
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FIG. 3: Radial wave function of a test particle falling radi- 
ally into a Schwarzschild black hole at -Robs = 7000M. We 
only show the I = m — 2 mode. The standard deviation pa- 
rameters a axe (a) l.OM and (b) O.lAf, respectively. The pa- 
rameters of oscillation frequency and the damping rate of the 
ringing tail are {Muj, M/r) = (a) (0.199, 0.015), (b) (0.314, 
0.063), respectively. 
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TABLE III: Comparison of our numerical results with the 
results of a free fall orbit {1 = 2 mode) 51]. 
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-Egw 
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relative error [%] 


5000M LOOM 


4.94 
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10" 


-3 
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5000 Af 0.50M 


8.38 
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10" 
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5000Af 0.25M 


8.78 
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10" 


-3 


4.57 


5000Af O.IOM 


8.96 
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10" 


-3 


2.61 


7000A4" LOOM 


4.17 


X 


10" 


-3 


54.7 


7000Ai" 0.50M 


7.94 


X 


10" 


-3 


13.7 


7000Af 0.25A4' 


8.71 
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10" 


-3 


5.33 


7000A'/ O.IOM 


8.95 


X 


10" 


-3 


2.72 



We compare our result with well-known one calculated 
in a frequency domain |5l| . The radiated energy E^^ j fo r 
the I = 2 mode is E^^ = 9.2 x 10-^{n^/M) in Ref. |5ll. 
Our results are summarized in Table IIIII We have the 
same behavior as in the case of a circular orbit that the 
error is small as the standard deviation parameter be- 
comes small. However, the relative error is larger than in 
the case of a circular orbit because a particle moves close 
to the black hole horizon, where the gravitational field is 
strong. 



We also investigate the convergence order of our code 
in the free fall orbit. First we fix a as 0.5M and 
varies Ar* between 0.1 A/ and 0.5M. The radiated en- 
ergy Eg^ at i?obs = 7000M are 9.65 x 10~^[^i'^/M] for 
Ar* = O.IM, 9.54 x 10~^[fi^/M] for An = 0.2M (rel- 
ative error from the radiated energy of Ar^ = O.IM is 
1.17%), 8.70x10-3[^VM] for Ar* = 0.4M (relative error 
9.92%), and 7.94 x 10-^[^^/M] for An = 0.5M (relative 



FIG. 4: Radiated energy of gravitational waves from a ra- 
dially in-falling particle into a Schwarzschild black hole for 
the I = 2 mode. Circle, square, diamond, and triangle denote 
a = l.OM, 0.5M, 0.25M, and O.IM, respectively. Dotted line 
represents the numerical result of Ref. |51| . 



error 17.7%). Therefore our code has second order con- 
vergence in space. Next we fix the relationship a = Ar* 
and varies a between O.IM and IM. Roughly speaking 
our code has a first order convergence in a, since the rel- 
ative error in Table ITTll decreases one half as a becomes 
one half. Note that the case a — l.OM may be out of the 
convergence region. 

We also set the observer quite far from the source in 
order to distinguish easily between the unphysical burst 
and the physical waveform. In fact, we integrate to com- 
pute the radiated energy in the range of 1600M < < 
3000M. Note that the radiated energy is irrespective of 
the variation of the range of integration. In Fig. ^ we 
also show the dependence of the I = 2 mode Eg^ on the 
location of the observer i?obs- We set the observer as 
Robs > 5000Af since the inappropriate initial condition 
disturbs the burst of the waveform. From Fig. ^ we find 
that the total radiated energy for a — l.OM and 0.5M 
are quite different from that found by Davis et al. [sij 
and this is irrespective of the location of the observer. 
The waveform with a large standard deviation parame- 
ter does not represent the quasi-normal ringing from a 
point particle. 

From the comparison of our two numerical results with 
the two cases, we confirm that our code has an ability to 
handle gravitational waves from a particle in a free fall 
orbit and in a bounded orbit. The convergence order 
of our code is first order in a and second order in Ar*. 
We also find that we should set the parameter sets as 
-Robs > 5000A/ and cr < 0.25M. 
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IV. GRAVITATIONAL WAVE FROM A DUST 
DISK 

A. Construction of a dust disk 

Here we explain our method for constructing a dust 
disk in Schwarzschild spacetime. We make five assump- 
tions for constructing a dust disk as follows: 

1. A disk has no self-interaction with each component; 
namely, a disk is composed of test particles. 

2. A disk is thin so that the only components are lo- 
cated in the equatorial plane of a Schwarzschild 
black hole. 

3. The density distribution for each radius of the disk 
is uniform so that it is considered as a uniform ring 
for each radius. 

4. Each radius of the ring has a Keplerian orbit, that 
is to say, a circular orbit. 

5. A disk is dynamically and radially stable. 

From Assumption |31 our disk is only located outwards of 
the innermost stable circular orbit of a test particle in 
Schwarzschild spacetime. We show a brief picture of our 
disk in Fig. [S] Note that the disk is composed of test 
particles (from Assumption^ in a circular orbit at each 
radius. Since we have a symmetry in the azimuthal direc- 
tion of the particle motion, we can construct a ring from 
test particles using Assumption 13 Then, we compose a 
disk from rings. 



We regard the mass distribution per unit width of the 
disk at radius r as m(r) given by 

m(r) = 27rCMr^-", (4.1) 

where C , ^, and n represent the normalization constant 
that has a dimension of r"~^, the rest mass of a test par- 
ticle, and the factor of the mass distribution, respectively. 
Note that n = denotes the constant mass density dis- 
tribution of the disk. The total mass of the disk /idisk 
is 

Mdisk ^ dr m{r) 
JeM 

^ f ^,7TC^i[R^- - (6M)2-"] (n ^ 2). 
\2^^C^l[lIi{R/M)-\Il6] (n = 2). 

In order to construct a disk numerically from a ring, 
we set each ring at the location which satisfies the energy 

H = El + j^{E^' - ^p) a = 1, • • • , ^d), (4.3) 

where is the energy at r = 6Af, E^'^ is the energy at 
r — R, R is the radius of the disk, and is the number 
of the rings. 

In order to initiate the radial motion of a dust disk, 
we slightly deplete the angular momentum of each ring. 
Since we maintain the energy conserved during the de- 
pletion, we still follow the geodesic equation for a ring 
after the depletion. 




FIG. 5: Illustration of our disk composed of test particles 
in a Schwarzschild black hole. EH and R represent the event 
horizon and the outermost radius of the disk, respectively. 



B. Gravitational wave from a dust disk 



Since we neglect the interaction between each compo- 
nent of the disk, gravitational waves emitted from the 
disk are obtained by superposing the waves emitted from 
each ring. For each gravitational wave emitted from each 
ring, we impose the same initial condition (Eq. 13.11 ) and 
the same boundary condition (Eqs. |3.2| and |3.3| ) as we 
used for the case of a test particle. In order to avoid the 
wave reflection at the boundary after a long evolution, 
we take the spatial grid as -5000Af < < lOOOOM. 
We also set the observer at i?obs — 5000Af . The radial 
wave function of gravitational waves \I>|'^"^'^-' is given by 

^f^^)^t, r) = ^ r dro m(ro) v&l^^'^/— ) 

Mdisk JeM r(0)=ro 

(4.4) 

where tq represents the radius of the ring at t — 0. Only 
the m — mode contributes to the gravitational waves 
from the ring due to its axisymmetric nature, where m is 
an index of spherical harmonics Yi^. Therefore the am- 
plitude of the gravitational wave is zero until the angular 
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t-r[M\ 
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t-r[M\ 



FIG. 6: Gravitational waveform from a dust disk at i?obs = 
5000M in the equatorial plane with constant density for four 
different disk radii, i.e., R = (a) IQM, (b) 14M, (c) 20M, and 
(d) 30A/. 



momentum is depleted. We only take the I = 2 mode 
into consideration because it is the dominant mode for 
the observer in the equatorial plane. In this paper we 
consider only radial mode due to m = 0. However the 
radial mode has a fundamental physics and there is also 
in a realistic situation. Thus it is considered that we can 
extract some important evidence by this analysis. 

We define the power spectrum of the gravitational 
waveform as 



P(lj) = lim — 



T/2 



T/2 



(4.5) 



where cj is a frequency. In our computation, we choose 
as large a T as possible, namely T = 4000M, assuming 
that the phase difference apart from the periodic waves 
at the edge of the integration should be negHgible. Note 
that we do not even integrate one long period of the wave 
in Fig. IHl(a) and Fig. 0(a) (b) (c). The valid frequency 
range of the power spectrum is several x (2tt/T) ^ w, 
namely 0.005 ^ Mu, in all our computations. Hereafter 
we calculate the gravitational waves with Ar* = a = 
0.25M and At = 0.5 x Ar*. 



1. Constant rest mass density distribution 

First we assume that the rest mass density distribu- 
tion of the disk is uniform, namely n = in Eq. H4.1|l . 
We vary the outermost radius of the disk as i? = lOM, 
14M, 20M, and 30Af . We deplete 1% of the angular 



momentum of each ring at {t — r) /M = to initiate the 
radial motion. As a matter of fact, the innermost radi- 
ally stable fragments of the disk approach r 7.85M. 
We show the gravitational waveform from a dust disk in 
Fig. |H| Since a uniform ring with a circular orbit retains 
an axisymmetric nature in the system, the system does 
not emit gravitational waves before the depletion. 
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FIG. 7: The power spectrum of the gravitational waveform in 
Fig.ini We extract the time range of lOOOM <t-r< 5000M 
to compute the power spectrum. The radius of the disk R is 
(a) lOAf, (b) 14A/, (c) 20M, (d) 30M. The values correspond 
to the peaks of the spectrum. 



We choose iVd = 50, 100, 175, and 300 for R = lOM, 
1AM, 20M, and 3QM, respectively, so that the step-size 
of the integration in Eq. (|4.4() should be the same. As the 
distance between the closest two radii of the components 
at t = is 0.035M ~ 0.26M depending on the radius, the 
components have sufficient numbers of particles to com- 
pose a disk. The first burst in the waveform is regarded 
as an inappropriate choice of the initial condition, as we 
mentioned in Sec. IIIII 

Since the waveforms (except for the disk with R = 
lOAf) are very similar to each other, it is difficult to ob- 
tain information on the disk solely from the waveforms. 
However, we find that there is a clear difference in the 
power spectrum of a gravitational waveform, varying the 
radius of the disk (Fig. [7|). Note that the spectrum of 
Fig. 13 is calculated from the waveforms of Fig. |2| To 
construct a power spectrum, we take the data set of a 
waveform in the range of lOOOM <t-r< 5000M. The 
value in Fig. corresponds to the frequency Muj of each 
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TABLE IV: Relationship between the initial radius of the ring 
and the characteristic frequency of gravitational waves from 
a ring. 

Rrins" [M] MuJi^'' 

8.0 0.0204 

10.0 0.0204 

14.0 0.0141 

20.0 0.00942 

30.0 0.00471 



'^iJring: Radius of the ring 

'^Mujgvj- Characteristic frequency of gravitational waves from a 
ring induced by 1% depletion of the angular momentum. 

peak. From these power spectra, the peak frequencies are 
almost the same, irrespective of the disk radius. How- 
ever, for R = 14M, 20Af, and 30M there is another peak 
frequency, which is lower than the frequency of the first 
peak in the case where R — lOM. 



0.025 




6 10 14 18 22 26 30 



Rring [M] 



0.0204 corresponds to the characteristic frequency of a 
particle at the outermost edge of the disk. Note that the 
peak frequencies above Mlo — 0.0204 correspond to the 
higher-order frequency of Mcu = 0.0204. In fact, the cor- 
responding frequencies have the following rule: 0.0408 = 
0.0204 X 2, 0.0628 w 0.0204 x 3, 0.0832 w 0.0204 x 4. 

The peak of the frequency below Muj — 0.0204 comes 
from the ring at the outer edge of the disk, and the peak 
at Mui = 0.0204 comes from the ring at the inner edge of 
the disk. The reason for no peak at the frequency below 
Muj = 0.0204 of a disk with R = lOM is that the fre- 
quency corresponding to R— lOM is the same as that of 
the peak at Moj = 0.0204. The phase cancellation effect 
plays a role in the emission of gravitational waves from 
the rings between the two edges of the disk. Therefore, 
it is possible to determine the radius of a dust disk from 
the power spectrum if the disk radius is ^ 12M. 

The characteristic frequency ujg„ coincides with the ra- 
dial frequency ujp of a particle ai r — R^ng due to the 
axisymmetric motion of the disk that corresponds to the 
TO = mode. In Fig. |S1 we plot the relation between 
the radial frequencies cUp of the particle at i?ring and ini- 
tial particle radius i?ring- Note that we deplete 1% of 
the angular momentum of the ring to initiate the radial 
motion. Since our frequency resolution is MAuj w 0.005 
as we mentioned before, the frequency cOp and LUg„ have 
a good correspondence. As mentioned above, because 
a peak corresponds to Wgw at i?ring = R in the power 
spectrum of gravitational waves, it is expected that we 
can determine the radius of a dust disk using the relation 
between ujp and i?ring- The quasi-normal ringing in the 
waveform represents the character of black hole; we can 
determine the character of the central black hole from the 
ringing. In a dust disk system, however, we cannot find 
the quasi-normal ringing in the waveform because there 
is a continuous inflow of fragments into a Schwarzschild 
black hole. Therefore the ringing is canceled by the differ- 
ent phase of the wave generated by each inflow fragment. 



FIG. 8: Radial frequency cup of the particle as a function of 
initial location i?ring. 



In order to understand the peak frequencies in the 
power spectra, we calculate the characteristic frequency 
of a particle around a Schwarzschild black hole. We 
summarize the characteristic frequencies of gravitational 
waves from a test particle with 1% depletion of the angu- 
lar momentum from a circular orbit in Table HVl From 
Fig. Hand Table Hvl the peak frequency Mlj = 0.0204, 
which appears in all disk models, corresponds to the char- 
acteristic frequency of a ring Wgw of i?ring = 8M and 
lOAf, and another peak frequency lower than Mlo = 



2. Non-constant rest mass density distribution 

Next we consider the non-constant mass density dis- 
tribution of a dust disk, which means that the rest mass 
density of the outermost edge of the disk is lower than 
that of the innermost edge, i.e., n > 0. In Fig. |51 we 
show a waveform from a dust disk with R = lOAf and 
R = 30M for n 0.5, 1.0, and 2.0, respectively. For 
a disk with R — lOM, the waveform does not strongly 
depend on the distribution. For the disk models with 
R — 30 Af, the waveforms approach those of the disk 
models with R = lOM as n becomes larger. Therefore it 
is also difficult to obtain information on the disk from the 
observational waveforms. However, in the power spectra 
(Fig. I10|l , the frequency corresponds to the peak having 
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t-r[M] t-r[M] t-r[M] 

FIG. 9: Gravitational waveform of a dust disk as a function of t — r. The parameters (n, R) are (a) (0.5, lOM), (b) (1.0, lOA/), 
(c) (2.0, 10M),(d) (0.5, 30M),(e)(1.0, 30M), (f) (2.0, 30M). We set the observer in the equatorial plane at iiobs = 5000Af. 
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FIG. 10: Power spectrum of gravitational waveforms in Fig. |2l We extract the time range of lOOOM t — r < 5000Af to 
compute the power spectrum. The values correspond to the peaks of the spectrum. 



the same feature as the one in the case where n = 0. 
Thus, the peak frequencies in the power spectrum cor- 
respond to those of gravitational waves from a ring at 
i?ring = R and at the inner region of the disk. Therefore, 
we conclude that we can determine the radius of the disk 
irrespective of the density distribution using the power 
spectrum of gravitational waves. 



3. Different type of angular momentum depletion 

Here we show our result based on a different type of 
angular momentum depletion of test particles whether 
our statement largely depends on the type of depletion. 
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FIG. 11: Radial frequency cup of the particle as a function of 
initial location i?ring- Dash line corresponds to Up of the par- 
ticle with a different angular momentum depletion (Eq. 14.61 ) 
and the solid line with the same angular momentum depletion 
as Fig. IHl 
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FIG. 12: Power spectrum of gravitational waves from a dust 
disk with a different type of angular momentum depletion 
(See Eq. 14.61 '). We choose the density profile of the disk as 
uniform, with the disk radius R = 20M. The values in the 
figure correspond to the peaks of spectrum. 



We deplete the angular momentum of the fragments in 
the following manner: 



The gravitational wave spectrum from the constant rest- 
mass density distribution is shown in Fig. 1121 Although 
the orbit of each fragments has been changed (the char- 
acteristic frequency is shown in Fig. 1111) . the peak fre- 
quencies in the power spectrum correspond to those of 
gravitational waves from a ring at i?ring = R and at the 
inner region of the disk. Therefore, we can determine the 
radius of the disk in a different type of angular momen- 
tum depletion using the power spectrum of gravitational 
waves. 



V. CONCLUSIONS 

We study gravitational waves from a dust disk around 
a Schwarzschild black hole using a black hole perturba- 
tion approach in a time domain. We especially focus on 
whether we can obtain the radius of the disk from grav- 
itational waves. 

We find that it is difficult to obtain information on the 
dust disk only by the observational waveforms. However, 
it is possible to determine the radius of the disk using 
the power spectrum of gravitational waves, irrespective 
of its density distribution. Our conclusion is quite simi- 
lar to the findings of Saijo and Nakamura 44, 45j. They 
claim that when a star is tidally disrupted by the central 
rotating black hole, we can determine the radius of the 
star from the spectrum of gravitational waves. In addi- 
tion, there is also a peak corresponding to a frequency 
of gravitational waves from a particle located at the in- 
ner edge of the disk. Since the standard deviation of the 
inner edge of the disk from R — 6M is a consequence 
of the radiated angular momentum, we could estimate 
the amount of radiation from the radius standard devi- 
ating from the innermost circular orbit. Although the 
detail disruption process requires 3D general relativistic 
calculation with the whole energy transport process, our 
statement from the simple model is still helpful to un- 
derstand the disruption process by observing the gravi- 
tational wave spectrum. 

We also mention the target of the gravitational wave 
source and its detectability. The typical frequency and 
the strength in our dust disk model are 
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Therefore it is possible to detect gravitational waves from 
a dust disk oscillating around a supermassive black hole 
by LISA. 

We have only investigated the outermost radius of the 
disk up to ~ 30M due to a limitation on computational 
time. However, our finding can be extrapolated into a 
more general astrophysical situation, such as the black 
hole formation phase of a supermassive star collapse. In 
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this phase, the collapse forms a supermassive black hole 
and a disk [s^ [s^ , and the fragments of the disk could 
fall into the supermassive black hole due to some dissi- 
pative mechanism or instabilities of the fluid that lead to 
a different configuration. The interaction between each 
material components of the disk is considerably small for 
soft equation of state such as a supermassive star, the 
radial oscillation of the disk may take a dominant role in 
exciting a peak in the power spectrum of gravitational 
waves. Therefore it could be one source generation sce- 
nario of gravitational waves useful for determining the 
size of the disk. 
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